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Mesoscopic fluctuations of two-port conductance and four-port resis- 
tance between Hall plateaus are studied within a realistic model for a two- 
dimensional electron gas in a perpendicular magnetic field and a smooth dis- 
ordered potential. The two-port conductance distribution P(g) is concave 
between g = and g = 1 and is nearly flat between g = 0.2 and g = 0.8. 
These characteristics are consistent with recent observations. The distribu- 
tion is found to be sharply peaked near the end-points g = and g = 1. The 
distribution functions for the three independent resistances in a four-port Hall 
bar geometry are, on the other hand, characterized by a central peak and a 
relatively large width. 



I. INTRODUCTION 

The dimensionless conductance g of small metallic disordered systems at very low tem- 
perature (where the phase- coherence length is larger than the linear size of the system), is 
characterized by its distribution P(g) which is nearly Gaussian with mean value propor- 
tional to the size of the system and variance of order unity. These properties, also termed 
as universal conductance fluctuations, are valid deep inside the metallic regime, for which 
g ^> 1. The question whether such (or different) fluctuations occur also in the critical regime 
where g is of the order of unity has been addressed in connection with the Anderson metal 
insulator transition prevailing in three-dimensional disordered systems [JTJ . In this case the 
critical region can be approached from the metallic side using scaling arguments and expan- 
sion in 1/g. Unfortunately, this is not possible in the quantum Hall regime, since insulating 
phases exist on both sides of the critical point. 
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In a recent experiment the conductance distribution function at the quantum Hall 
transition point was elucidated. The geometry of the experiment was that of a two- 
dimensional strip with two-ports, namely, one source and one drain. The transition point is 
defined as that value of the magnetic field for which the averaged conductance is (n + 1/2)^-. 
(Since the focus of study is mesoscopic fluctuations, the notion of transition point is more 
appropriate than that of critical point, which, strictly speaking, applies only in the ther- 
modynamic limit). It was found that in the range of values 0.1 < g < 0.9 the distribution 
P{g) is roughly flat. Subsequent theoretical work was mostly based on the network model 
||. The qualitative features of the experimental distribution have been obtained by nu- 
merical calculations Q based on the network model with free boundary conditions on the 
transverse coordinate. Calculations were performed also with periodic boundary conditions 
in order to elucidate the role of edge states. Roughly speaking, adopting free (periodic) 
boundary conditions corresponds to calculating R xy (R xx ). Evaluation of P(g) using real 
space renormalization group methods applied to the network model resulted in a concave 
shape with a very flat bottom |5]|| and sharp peaks near g = and g = 1. It is therefore 
interesting to check (against other models) whether the network model which proves to be 
extremely successful in describing the critical behavior of the integer quantum Hall effect is 
also capable of encoding the correct statistical properties of the conductance for mesoscopic 
systems. 

A natural extension of the physics discussed above is related to the mesoscopic fluctua- 
tions of the Hall resistance. To the best of our knowledge, investigation of the Hall resistance 
fluctuations in the transition points between Hall plateaus has not yet been reported. It 
might be worth pointing out here that the precise relation between the two-port conductance 
g, discussed in connection with the experiment of Ref. 0, and the conductivity tensor a^, 
defined through the Kubo formula, is not unambiguously clear. In particular, an interesting 
sum-rule relates the two-port transmission coefficient to a combination of Hall and dissipa- 
tive resistivities J7|. This relation holds, in fact, not only at the transition region between 
plateaus. In order to separate the effects of the Hall and dissipative resistivities a multi-probe 
measurement is required, as employed in the original experiment ||. Once a multi-probe 
geometry is defined, numerous resistances can be determined (beside the Hall resistance). 

In order to address these questions, calculations of conductance (resistance) distributions 
between Hall plateaus for two (four) port setups are presented below (in the four-probe 
geometry, the appropriate quantity is in fact the resistance). The Schrodinger equation 
for an electron moving in part of the plane restricted by certain boundaries and subject 
to a strong perpendicular magnetic field and a smooth disordered potential is solved, and 
transmission and reflection coefficients between various ports are computed PJT0|] . The two- 
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port conductance is calculated using the Landauer [|ll]] formula, while the various four-probe 
resistances are calculated from the Biittiker formula ||12|| . Note that no projection on lowest 
Landau level is made. 

We find that the distribution P(g) of the two-port conductance is concave, with rather 
shallow minimum at g — 1/2. As we have already pointed out, the experimental data in- 
dicate a flat distribution between g = 0.1 and g = 0.9, and in this domain the calculated 
distribution is quite smooth and hence qualitatively accounts for the experimental distribu- 
tion. On the other hand, there is so far no experimental evidence for the occurrence of two 
peaks near g = and g = 1. It is also found that conductances g{Ep, H) at different points 
of the fermi energy (E) and magnetic field (H) are correlated if the two points are along a 
transition line. 

In the four-probe geometry there are three independent resistances that can be measured, 
i?i2,34 -Ri3,24 and -Ri3 5 i3. In this notation Rij^i is the ratio of the voltage difference between 
ports i and j and the current flowing between ports k and I. Inspecting the histograms of 
these three independent resistances at the transition point (defined by < -Ri3,24 >= 0.6-^) > 
show that all three of them have a distribution which is, more or less, close to a normal 
one. The width of -Ri2,34 (the so called bend resistance) is of order AR = 0.2h/e 2 while 
for -Ri3 5 24 (the Hall resistance), and i?i 3> i 3 the order of the width is AR = 0.5/i/e 2 . This is 
quite large compared with the corresponding average. Thus, within the realm of mesoscopic 
fluctuations, the Hall resistance is not a universal number. We are unaware of an experiment 
aimed at studying mesoscopic fluctuations of multi-probe resistances. 

The method of calculation is briefly explained in Sec. II, and the results are presented 
and discussed in Sec. 



II. FORMALISM 

The two-port system is schematically displayed in Fig. |l]a. It consists of a strip — 00 < x < 
00, —L/2 < y < L/2 which confines the motion of electrons whose wave function vanishes 
on the walls. The four-port system is schematically displayed in Fig. [l|b. It consists of 
two perpendicular strips of width L which intersect each other. It is convenient to number 
the ports (leads) from the left in a counter clockwise sense, so lead number 1 is the left 
lead etc.. In both two and four-port systems, the motion of the electron is ballistic except 
within the domain (—L/2 < x < L/2, —L/2 < y < L/2) where the electron is subject to a 
smooth disordered potential, defined below. The system is subject to a strong perpendicular 
magnetic field. It is useful in this geometry to define the vector potential A in the Landau 
gauge, A x = —Hy, where H is the magnetic field strength On this square one considers a 



3 



lattice of N 2 Gaussian potentials of width a, centered at the points (x n ,y m ) = ( — L/2 + 
a/2 + (n — l)a, —L/2 + a/2 + (m — l)a), with n,m = —N/2 + 1, • • • , N/2, assuming that 
iV = L/a is an even integer. The random potential is written as 

V( X ,y) = ^ nmeXp (- <- X - X ^ 9 + ^-^ ). (1) 

nm \ ^ J 

The impurity strengths v nm are independent random energies, distributed (for convenience) 
uniformly between —w/2 and w/2. For this distribution one has 

(Vnm) = 0, (2) 

and 

w 2 

{^nm^n'm') "^^nrn'^mm' • (3) 

The single-particle reduced Hamiltonian then reads 
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+ V{x,y)=H + V. (4) 



The magnetic length is given by £h = (hc/eH) 1 ^ 2 , where — e < is the charge of the electron. 

In order to calculate transmission and reflection coefficients one needs to obtain the 
asymptotic states away from the central region, i.e., in the leads. Consider for example the 
left lead — oo < x < —L/2, —L/2 < y < L/2. For a given Fermi-energy Ep the solution of 
the free Schrodinger equation Tio^ = Ep^ in this domain is given by 

^(x,y) = q-^e^+^fiy; ±q n ), (5) 



where f(y; q) is a linear combination of parabolic cylinder functions [0 such that the 
boundary conditions f\—L/2;q) = f(L/2;q) = determine the eigen-momenta q n {Ep,H). 
Note that in the case of a four-port system, the role of x and y is interchanged in the vertical 
leads. For convenience, the amplitude of the incident plane wave is set to be exactly unity 
on the boundaries between the ideal conductors and the region of impurities. 

It is worth mentioning that for the present geometry, it is not possible to impose periodic 
boundary conditions in the transverse direction. This situation is distinct from the one 
encountered in the network model ||. 

The number M of open channels is the maximum number of positive momenta. It is an 
increasing number of the Fermi-energy Ep and a decreasing number of the magnetic field 
H. The behavior of the momenta q n (Ep,H) as function of the two parameters has been 
discussed elsewhere M. 
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Consider now a solution of the full scattering problem such that there is an incoming 
wave ipt °f un ^ amplitude in one of the leads (say lead i), and outgoing waves in all leads 
j (including the lead j =i). The wave function in lead j is then given by, 



M 



K=J2^n,(forj^i), (6) 



m=l 

and 

M 



\£? u = ih+ + y ih~t u . (7) 



m=l 



The matrices t^ n and are transmission reflection amplitudes respectively. Within linear 
response formalism, the two-port conductance g and the four-port resistance tensor Rij^i can 
be determined by these amplitudes. In the first case the relation is given by g = Tr[t 12 t 12 ^} 
while in the second case the relation is given in Ref. \I2\. The essential part is: Rij^i oc 
(T kl T lj - T kj T li ), where T ki = Tr[t ki t k % 

We recently introduced a variational algorithm for the generalized Z-matrix which, given 
a finite region in configuration space, effects the following transformation of the scattering 
wave-function on the boundary surface: 

Z{aV{S) + bDV(S)) = dV{S) + qDV(s). (8) 

In the above we define the generalized normal derivative D^(S) = n • ^tt^(S), where 
7r = p — is the generalized momentum for a particle of momentum p and charge e in 
a vector potential A. The parameters a, b, d, and q are arbitrary except that they obey 
the constraint aq — bd = —1. In all studies we use a = —q = 1, b = d = 0, for which Z 
is the negative of the log-derivative matrix. Actually, we use the inverse of this matrix, as 
discussed in reference WM , which is the negative of the well- known R matrix. The notation 



^(S) indicates that the wave-function is evaluated on the boundary surface S and projected 
onto a suitable basis defined on the surface. In all studies we use sine functions on the 



boundaries of the regions in Fig. 1. Using the coefficients of f(y; ±g n ) in the sine basis, [TJ 
we can plug Eqs. (6) and (7) evaluated on the boundaries of the scattering structures into 
Equation (|S|) and solve for the coefficients using standard singular-value decomposition 
algorithms |16| . 

Because of the large bases required, direct computation of the R-matrix was not conve- 
nient. Instead, we used a straightforward extension of a recursive cellular R-matrix algorithm 
outlined by Nesbet [ 17f . In the 2-port case, we divided the coordinate x into N x subregions. 



On each subregion, we computed an R-matrix using a basis consisting of sine functions in 
y which go to zero at y = ±~ and a Legendre discrete variable representation (DVR) in 
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x Jl8| (for a discussion of the Legendre DVR we used, see Ref. |I4|). Using the continuity 
and differentiability requirements on \1/ at the boundaries between cells, we concatenated 
adjacent R-matrices into global R-matrices for the entire region, as discussed by Nesbet. 

In the 4-port case, we divided the region in Fig. lb into N% square subregions; that is, we 
divided x and y each into N s subregions. On each subregion we computed an R-matrix as 
discussed in Ref. [II] and concatenated them into the global R-matrix for the entire region 
in a procedure analogous to that outlined above. A few cautions may be in order here. 
First, for the subregion R-matrices we use Legendre functions as surface functions, not sine 
functions. This modification is necessary since to use sine functions on the subregions would 
be to imply that the scattering wave-functions are zero at all the points of intersection 
between any four cells. Second, the scattering wave-function is zero at the four corners 
of the square in Fig. lb. These conditions are enforced by means of linear constraints on 
the wave-functions in the subcells at these corners. Finally, since the final global R-matrix 
computed in the manner outlined above would be expressed in terms of Legendre functions 
confined to 4N S subregions on the surface, the resulting matrix must be projected onto the 
desired surface sine functions. 



III. RESULTS 

Let us now make a choice of the potential parameters. We first assign a relation between 
a and a such that the exponential factor in the overlap integral between Gaussians of two 
adjacent impurities, calculated to be 27rcx 2 exp (— ^sr) equals exp(— 1/2). Hence, a = a/y/2. 
Next, the length of the system L and the density of impurities p = (Na/L) 2 are determined 
by the requirement that in the range of pertinent magnetic fields, the length of the system is 
much larger than the magnetic length and there is about one impurity per magnetic area i\ . 
We have found that computation time grows quickly with the total number of impurities, 
N 2 . The present calculations are carried out with iV = 10. This yields L « 10£#. 

Hereafter we use atomic units; length is given in units of the Bohr radius r = 0.5292 
Angstrom, energy is expressed in Hartree (1 Hartree (e 2 /r ) = 27.2 eV), and the atomic 
unit of magnetic field is equivalent to 1715.1 Tesla. Calculations are performed for system 
width L = 2078 Bohr. The value of disorder strength (w) must be commensurate with 
those of Fermi-energy and magnetic field so that the conductance is mainly dominated 
by magnetic field effects and not by disorder. We chose w = 4.26 x 10~ 5 Hartree and 
require that the Fermi-energy Ep and the magnetic field H will vary in the region for 
which the transition between the first plateaus to the insulating phase occur. These are 
1.3 x 10~ 5 < E F < 1.5 x 10~ 5 and 3.45 x 10~ 3 < B < 3.80 x 10~ 3 respectively. 
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A. Two-Port System 



First we show in Fig. |2] the two-port average conductance together with its variance 
near the last plateau. This figure shows that mesoscopic fluctuations occur mainly in the 
region of transition between plateaus, in agreement with experiment ||. The distribution 
P{g) of the two-port conductance is displayed in Fig. |3| for three values of the magnetic 
field for which the average value of the conductance is 0.514, 0.498 and 0.478. In all three 
cases the distribution is concave with very shallow minimum at the center and maxima at 
the endpoints g = and g — 1. The region between g = 0.2 and g = 0.85 is roughly 
flat, in agreement with experiment M. Note that the distribution presented here is similar 
to the one predicted from the network model using the renormalization group technique 
P||. On the other hand, its characteristics near the end-points is distinct from the one 
obtained by numerical solution of the network model (see Fig. 5b in Ref. [|J). Assuming the 
latter procedure to be more accurate than the renormalization group one, we suspect that 
the network model and the present strip geometry model yield distributions which are very 
similar in the main part of the interval < g < 1, but are disparate around the end-points. 
On the other hand, the results of the present model are in good agreement with RG analysis 
of the network model. The question which one is closer to the real distribution requires 
further experimental study. 

Another question which is interesting in this context is whether conductances g(Ep, H) 
at different points in the (Ep, H) plane are correlated. Preliminary results indicate that 
the answer to this question is affirmative p0| . A useful concept in mesoscopic physics is 
that of ergodicity, which roughly speaking, implies that fluctuations (of various quantities 
such as the conductance) as function of a given parameter are equivalent to sample-to- 
sample fluctuations. The experiment discussed here M is based on this hypothesis, since 
measurements are performed on a single sample at many values of the magnetic field. In other 
words, it is tacitly assumed that, in the thermodynamic limit, (g(Ep, H)g(E' F , H') decreases 
very rapidly with increasing \Ep — E' F \. The question that can be posed here is whether 
nontrivial correlations (g(Ep, H)g(E' F , H')) do exist. Within the limitation on system size 
as explained above it is not possible to give a definite answer here. Instead, we display in 
Fig. |] a contour plot of dg/dH in the (Ep, H) plane. The regions of maximal derivatives are 
clearly shown in the projected contour plot. In the bulk geometry they simply correspond to 
the Landau levels E n = (n + \)Tiui. Here they correspond to transition between n and n + 1 
edge states. It might then be argued that if (E 1 ,H 1 ) and (E2,H 2 ) are two different points 
along one of these lines then the correlation (g(Ei, Hi)g(E 2 , H 2 )) is larger than otherwise. 
Yet, this analysis is only qualitative. In particular, the experimental analysis of conductance 
derivative reveals correlation lines which cross each other ||20|| . Interpretation of these latter 
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observations apparently goes beyond the single particle picture. 

B. Four-Port System 

We now present results for resistance distributions in the four-port geometry. As men- 
tioned in the introduction, there are three independent resistances to be studied, #13,24 (the 
Hall resistance), #12,34 (the bend resistance) and #13,13. Note that the latter should not be 
confused with the dissipative conductance R xx since current and voltage are measured at the 
same points in the leads. Calculations (and measurement) of R xx require a six port geome- 
try. This is the reason why, in the present scheme, we can calculate the Hall resistance and 
not Hall conductance. Nevertheless, the resistance #13,13 has a non-trivial physical content 
of its own |T0| . 

The value of the magnetic field H c at which the quantum Hall transition occurs is usually 
determined by that value for which R xx has a sharp peak, or, equivalently, the Hall conduc- 
tance takes on the values G xy = n + 1/2. The corresponding value of the Hall resistance 
Rxy = r2 l v r 2 at the transition point depends on the corresponding value of G xx . Assuming 
G XX (H C ) = 0.5 ||19|| , implies that at the transition point between n = 2 (where R xy = ^2) 
and n = 1 (where R xy = ^) one has R xy (H c ) = 0.6^ (and not 0.75^- as one might be 
tempted to think). We choose to present our results at this point and not the last transition 
point since there R xy = ^2 exactly as in the plateau. 

In Fig. ^ we display the distributions of the three resistances at H c . All the resistances 
are characterized by a central peak and small tails. In all three cases, the width of the 
distribution is of the same order as the average itself, which means that mesoscopic fluctua- 
tions are rather significant. We are tempted to conjecture that the distribution of the Hall 
conductance would be similar to that of the Hall resistance, which, comparing with Fig. |3] is 
much distinct from that found for a two-port system. This is quite significant since, based 
on the edge state picture, it is intuitively suggestive that two-port magneto-conductance at 
strong magnetic field is the same as the Hall conductance, which just counts the number 
of edge states. This is definitely true in the plateau regions, and apparently also in the 
transition region if the average quantity is considered. On the other hand, the pertinent 
mesoscopic fluctuations are distinct. We have yet to find an explanation for this difference. 
Unfortunately, measurements of resistance fluctuations in multi-port systems are not yet 
available. 

In conclusion, we have presented a study of conductance and resistance distribution 
functions in the quantum Hall regime, using a realistic calculation scheme. The two-port 
conductance distribution agrees well with the experimental one 0. Its behavior is charac- 
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terized by a concave shaped curve which is rather flat at the center and sharply peaked near 
the end-points g — and g — 1. Our result also indicate possible correlations of conduc- 
tances at different points in the (E, H) plane, if these points are located on a transition line. 
The four-port resistance distributions are peaked at their respective central values and have 
large widths (compared with their mean values). It remains to be determined whether these 
fluctuations survive in the thermodynamic limit at zero temperature. 
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IV. FIGURE CAPTIONS 

FIG. 1. Sample geometry used in the present calculations, a) Two port system and b) four-port 
system. The shaded area contains a disorder potential. 

FIG. 2. two-port average conductance (solid line) and its variance (dash-dotted line) as a func- 
tion of the magnetic field for fixed Fermi energy E = 3.65 x 10~ 3 Hartree. The system width is 
L = 2078 Bohr, and the disorder strength is w = 2.3 x 10 -5 Hartree. 
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FIG. 3. Distribution P(g) of two-port conductance in the middle of the transition between Hall 
plateaus n = 1 and n = for three different values of the magnetic field, B\ = 0.00360 a.u., 
B<i = 0.00365 a.u. and B% = 0.00370 a.u. The corresponding values of the average conductances 
are g± = 0.514, gi = 0.498 and 53 = 0.478. Parameters are as in figure ^. 

FIG. 4. Three dimensional and contour plot of the conductance derivative dg(Ep, H)/dH in 
the (Ep,H) plane for a given sample. Length and disorder parameters are as in Fig. g. Units of 
energy and magnetic fields are specified in the text, while g is in units of e 2 /h. 

FIG. 5. Displayed from top to bottom are distributions of Hall resistance -^13,24, bend resistance 
-^12,34 and the "through" resistance -^13,13 for H = 0.00135 a.u. for which the average value of the 
Hall reistance is 0.6. 
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